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We explore the thermodynamics and phase structure of the Polyakov loop-extended two flavor 
chiral quark-meson (PQM) model beyond the mean-field approximation at imaginary chemical 
potential. Our approach is based on the functional renormalization group (FRG) method. At finite 
temperature and imaginary chemical potential, we solve the renormalization group flow equation 
for a scale-dependent thermodynamic potential in the presence of the gluonic background field. We 
determine behavior of order parameters of the PQM model in the FRG approach and compute the 
phase diagram. We compare our FRG results with that obtained in the mean-field approximation 
at imaginary chemical potential. 
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I. INTRODUCTION 

Thermodynamic properties of strongly interacting 
matter at nonzero baryon density and at finite temper- 
ature have been explored numerically within the Lattice 
Quantum Chromodynamics (LQCD) The LQCD 

results show that QCD at physical quark masses and at 
vanishing baryon density exhibits restoration of the chi- 
ral symmetry and deconfinement at T ~ 160 MeV 
Q . The crossover nature of the transition at zero baryon 
density makes the variation of the order parameters quite 
smooth and leads to the broader width of the pscudocrit- 
ical temperature Q- 

Unfortunately, the thermodynamics of strongly in- 
teracting matter at large baryon densities is presently 
not accessible in the first principle LQCD calculations, 
because of a complex structure of fcrmion determi- 
nant. Several methods circumventing this, so-called, 
sign-problem have been introduced. However, they are 
limited to small values of chemical potentials or small lat- 
tice sizes [l^l- Thus, LQCD calculations have not been 
able yet to address a nature of a critical behavior at finite 
density on sufficiently large lattices and at small quark 
masses. 

Phenomenological models and effective theories offer a 
viable framework for exploratory studies. The properties 
of low-energy hadrons as well as the nature of the chiral 
phase transition at finite temperature and density have 
been studied intensively in such effective models [ll| - l30| . 
Recently, the physics of color confinement and its relation 
to the chiral symmetry breaking have been also addressed 
in such a framework. Particularly successful, are models 
obtained by extending the chiral Lagrangians, such as the 
Nambu-Jona~Lasinio or the quark-meson, by introduc- 



ing a coupling of quarks to uniform temporal background 
gauge fields, the Polyakov loop [l5|, [1^. Though, these 
models have been formulated so as to incorporate essen- 
tial features of QCD, nevertheless, properties of phase 
transitions at large baryon densities strongly depend on 
their parameterizations [H, [1^, [Sll - lssj . 

The imaginary chemical potential formulation pro- 
vides complementary constraints on the phenomenolog- 
ical models. In addition, at imaginary /x, the phase di- 
agram can be directly computed in LQCD by virtue of 
the absence of the sign problem. Within LQCD, these re- 
sults have been used to extrapolate the phase boun dary 
from an imaginary to a real chemical potential |34| - [39| . 
Such a strategy was found to be successful in two-color 
QCD [i^l, resummed perturbation theory [4l|, and quasi- 
particle models [42j . 

Another interesting issue, that can be addressed at 
imaginary chemical potential, is the order of the tran- 
sition at Roberge- Weiss (RW) endpoint, which might 
dictate the property of the transition at vanishing and 
finite real chemical potential. Recently, it was shown, 
within lattice simulations, that the critical behavior at 
the RW endpoint has a non-trivial dependence on the 
quark masses for both, two-fisj and three-flavor QCD 
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The formulation of the theory at imaginary /i also 
opens the possibility to obtain the canonical partition 
function from LQCD by means of the Fourier transfor- 
mation of the grand canonical partition function (45l - 
IZsj . Therefore, studies of thermodynamics at imaginary 
chemical potential are important to explore thermody- 
namics of the QCD at finite density and its phase dia- 
gram [igf . 

The Polyakov loop extended Nambu-Jona-Lasinio 
(PNJL) [itIi and quark-meson (PQM) ^ models re- 
produce essential features of the QCD thermodynamics 
already in the mean-field approximation at real and at 
imaginary chemical potential [13, HII- Especially, the 
Roberge- Weiss (RW) periodicity Z{fii/T) = Z{^i/T + 
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27rfc/iVc) and the associated Z{Nc) transition at high 
temperature [13 can be described in these models. 

In our previous studies [5l|, we have discussed the 
phase structure of the PNJL model at imaginary chem- 
ical potential. While the statistical confinement feature 
of this model naturally provides characteristic properties 
of the order parameters, the chiral condensate and the 
Polyakov loop, its thermodynamics, however, is strongly 
influenced by the choice of the Polyakov loop potential. 
Furthermore, to reproduce within the mean field approx- 
imation the LQCD results in terms of such effective mod- 
els, several amendments seem to be necessary [28l. [53l - l55| . 

To correctly account for the critical behavior and scal- 
ing properties near the chiral phase transition, it is neces- 
sary to go beyond the mean-field approximation and in- 
clude fluctuations and non-perturbative dynamics. This 
can be achieved e.g. by using methods based on the func- 
tional rcnormalization group (FRG) (56l - [6l| . 

In the following, we consider the PQM model at imag- 
inary chemical potential and study its critical properties 
and phase diagram within FRG approach. We formu- 
late and solve the suitably truncated FRG flow equation 
for fluctuations of the meson fields in the presence of 
the Polyakov loop which is treated as a background field 
on the mean-field level. We extend previous studies [HH 
to imaginary chemical potential and explore influence of 
fluctuations on the chiral and deconfinement order pa- 
rameters. We compare our FRG results at imaginary 
chemical potential with that obtained in the mean field 
approximation. We show, that there is an essential mod- 
ification of thermodynamics at imaginary chemical po- 
tential owing to quantum mesonic fluctuations. 

The FRG approach was previously applied to study the 
phase structure of two-flavor QCD in the chiral limit at 
imaginary chemical potential [49[. In the current paper, 
we perform calculations within the PQM model to study 
the role of the mesonic fluctuations on the structure of the 
phase diagram at imaginary chemical potential by com- 
paring results obtained within the FRG approach and 
undeer the mean-field approximation. 

In the next section, we introduce the PQM model and 
the implementation of the FRG method at imaginary ^. 
We present our results in Sec. IIIII Section |IV] is devoted 
to the summary. 



II. THE POLYAKOV-QUARK-MESON MODEL 

The quark-meson model is an effective realization of 
the low-energy sector of QCD, which incorporates chiral 
symmetry. Because the local color SUc{N) invariance of 
QCD is replaced by a global symmetry, the model does 
not describe confinement. Nevertheless, by introducing a 
coupling of the quarks to a uniform temporal color gauge 
field, represented by the Polyakov loop, some confine- 
ment properties can be effectively included [HI, [l^, [2^ . 



The Lagrangian of the PQM model reads [23| 



(1) 



The coupling between the effective gluon field and quarks 
is implemented through the covariant derivative 



(2) 



where Ai, — gA^\'^/2. The spatial components of the 
gluon field are neglected, i.e. Ai, = Si^qAq. Moreover, 
$*) is the effective potential for the gluon field ex- 
pressed in terms of the thermal expectation values of the 
color trace of the Polyakov loop and its conjugate 

<i> = i-(Tr,L(^)), <i>*^i-(Tr,it(^)), (3) 

with 



L{x) = ■Pexp 



I I drAi (x, r) 




(4) 



where V stands for the path ordering, (3 = 1/T and 
Ai — iAq. In the 0(4) representation, the meson 
field is introduced as = (c, 7?) and the correspond- 
ing SU{2)l ® SU{2)ii chiral representation is defined by 

a + IT ■ TT 75 . 

The purely mesonic potential of the model U{a,T:), is 
defined as 

Uia,Tr)^^{a^ +7f^ -v^Y -ca, (5) 

while the effective potential of the gluon field is 
parametrized in such a way as to preserve the Z{3) in- 
variance, 

(6) 

The parameters, 

b2{T) = ao + ar(|)+a.(f)%a3(|)' (7) 

with ao = 6.75, ai = -1.95, 02 = 2.625, 03 = -7.44, 
63 = 0.75, 64 = 7.5 and Tq = 270 MeV were chosen to 
reproduce the equation of state of the pure SUc{S) lattice 
gauge theory. When the coupling to the quark degrees 
of freedom are neglected, the potential ^ yields a first- 
order deconfinement phase transition at To- 

Alternative parametrization of the Polyakov loop po- 
tential, see e.g. Ref. [S^, niay provide a better fit of the 
lattice results by the PQM model. However, we follow 
our previous studies of the PQM model formulated at the 
real chemical potential, and apply the potential (H]) in the 
model calculations at imaginary ^. The main conclusions 
on the influence of the mesonic fluctuations on the phase 
structure will not change if different parametrization of 
the potential is considered. 
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A. The FRG method in the PQM model 

In order to account for mcsonic fluctuations in the 
PQM model, we employ a scheme based on the functional 
rcnormalization group (FRG). This scheme involves an 
infrared rcgularization of the fluctuations at a sliding 



momentum scale k, resulting in a scale-dependent efl^ec- 
tive action Ffe, the so-called effective average action [sgI - 
[5§|. We treat the Polyakov loop as a background field, 
which is introduced self-consistently on the mean-field 
level while fiuctuations of the quark and meson fields are 
accounted for by solving the FRG flow equations. 



Following our previous work [63 | , we formulate the flow equation for the scale-dependent grand canonical potential 
density fi^ = TTf^/V for the quark and meson subsystems at flnite temperature T and imaginary chemical potential 
/i = iOT as follows 



9fcr!fc($,$*;T, 



fc4 



127r2 1 



1 



l + 2nB{E„-T) 



(8) 



1 - 7V($, $*; T, e) - iV($, $*; T, ( 



Here nB{ET^^^;T) is the bosonic distribution function 

nBiE.^,a;T) 
with the pion and sigma energies 

E^ = ^Jk^+Ul , E„ = ^Jk^+T^l + 2pni 

where the primes denote derivatives oi CI — Q + ccr with respect to p field, p = {a^ + 7T'^)/2, and /3 = l/T. The fcrmion 
distribution functions iV($, T, 0) and 7V($, $*; T, 9), 



1 + 2$* e^pl^Eg - i0] +^e^p[2{l3Eg 



^ ' ' ' ' l + 3<S?cxp[2{PEg-ie)]+3'^*exp[l3Eg-i0]+exp[3{PEq-ie)]' 
iV($,$*;r,6l) iV($*,$;r, -6*), 

are modified because of the coupling to the gluon field. Finally, the quark energy is given by 



(9) 
(10) 



Eg - Vfc^+VP- 



(11) 



The minimum of the thermodynamic potential is de- 
termined by the stationarity condition 



da 



da 



-c = 0. 



(12) 



The flow equation (|8]) is solved numerically with the ini- 
tial cutoff A = 1.2 GcV (see details in Ref. The ini- 
tial conditions for the flow are chosen to reproduce the 
in-vacuum properties: the physical pion mass ttIt^ — 138 
MeV, the pion decay constant /^r = 93 MeV, the sigma 
mass rricr = 600 MeV, and the constituent quark mass 
niq = 300 MeV at the scale k ^ 0. The symmetry 
breaking term, c = tti^/tt, corresponds to an external 
held and consequently does not flow. In this work, we 
neglect the flow of the Yukawa coupling g, which is not 
expected to be significant for the present studies (see e.g. 

Rcfs. [ei,!!^). 



By solving the equation ([5]) , one obtains the thermody- 
namic potential for the quark and mesonic subsystems, 
ilk^oi^,^*;T,0), as a function of the Polyakov loop 
variables $ and $*. The full thermodynamic potential 
51($, $*; T, 9) in the PQM model, including quark, meson 
and gluon degrees of freedom, is obtained by adding the 
effective gluon potential W($, $*) to f2fe^o(^', 6): 



(13) 



At a given temperature and chemical potential, the 
Polyakov loop variables, $ and are then determined 
by the stationarity conditions: 



d_ 



17($,$*;T,e) = 0, 
■r2($,$*;r,6l) = 0. 



(14) 
(15) 
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The thermodynamic potential (|13p does not contain tion of the high momentum states by approximating it 

contributions of thermal modes with momenta larger as quarks interacting only with Polyakov loop degrees of 

than the cutoff A. We take into account the contribu- freedom as done in Ref. |62l |. 

I 



B. The mean- field approximation 

To show the importance of mesonic fluctuations on the thermodynamics of the PQM model formulated at imaginary 
fi wc compare the FRG results with those obtained in the mean-field approximation. In the latter, mesonic fluctuations 
are neglected and the mesonic fields are replaced by their classical expectation values. 

The thermodynamical potential of the PQM model in the mean- field approximation reads (23j . 

Qmf - Ui^, $*) + C/((a), (tt) = 0) + $*). (16) 

Here, the contribution of quarks with the dynamical mass iriq ~ g{a) is given by 

0,,((a),<i>,$*) = -^m^ln(^)-2iV^r J ^ {lng(+) ((a), 3>, $*; T, 0) + ln5(-)(('^), ^)} , (17) 

where 

5(+)((cr),$,$*;r,6l) = l + 3<i>cxp[-~{PEg-ie)]+3<i>*cicp[-2{PEg-ie)]+cxj>[-3{PEg^ie)l (18) 

g("'((CT),$,$*;r,6l) = .g<+)((cr),$*,$;r, -6*) (19) 



and Eq = y + "^g the quark quasi-particle energy. 

The first term in Eq. (|17p is a vacuum contribution reg- 
ularized by dimensional regularization with renormaliza- 
tion scale M [g^l- The relevance of the vacuum con- 
tribution for the thermodynamics of chiral models was 
demonstrated and studied in detail in Refs. [g^] and [68| . 

The equations of motion for the mean fields are ob- 
tained by requiring that the thermodynamic potential is 
stationary with respect to changes of (a), $ and <J>*. Uti- 
lizing the fact that is the complex conjugate of $ at 
imaginary chemical potential, we introduce the modulus 
and the phase of the Polyakov loop as $ = l^le'"^ and 
$* = l^je"**^ and the stationary condition is given as 

The model parameters are fixed to reproduce the same 
vacuum physics as in the FRG calculation. 

III. THERMODYNAMICS OF THE PQM 
MODEL 

The PQM model, which is expected to belong to the 
same universality class as QCD, exhibits a generic phase 
diagram with a critical point at a non-vanishing real 
chemical potential f23| . A detailed comparison between 
thermodynamic properties of the PQM model obtained 
within the FRG and in the mean-field approach at real 
chemical potential was recently studied in Ref. f69| . 

In the following, we explore the critical properties of 
the PQM model at imaginary chemical potential fi/T = 



—iO using the functional renormalization group approach 
and in the mean-field approximation. 

Under the mean-field dynamics, many features of the 
PQM model at imaginary chemical potential should be 
common with that obtained previously in the PNJL 
model. Thus, the nature of the phase structure in the 
mean-field approximation for the PNJL model found in 
Ref. [sif also applies to our present calculations within 
the PQM model. Therefore, we focus on the effects of 
the mesonic fluctuations on thermodynamics at imagi- 
nary chemical potential and at finite temperature. 

At imaginary chemical potential, in addition to the 
chiral transition which is characterized by the chiral or- 
der parameter a and the dcconfinement transition which 
is indicated by the rapid increase of the modulus of 
the Polyakov loop, the PQM model exhibits also the 
Roberge- Weiss (RW) transition. Since the thermody- 
namic potential of the PQM model shares the Roberge- 
Weiss periodicity of QCD, n{T, 9) = n{T, 6 + 2|ii) with 
n € N, which is a remnant of global Z{3) symmetry of 
the SUc(3) gauge group, the first-order RW phase transi- 
tion takes place in the deconfined phase at 9 ~ 7r/3. The 
RW transition can be understood as a sudden change of 
the phase of the Polyakov loop. While as a func- 
tion of 9 smoothly changes below the RW cndpoint, at 
T = Trw it has a discontinuity caused by a transition 
from one Z(3) sector to another. This transition also in- 
fluences the thermodynamic potential and the modulus 
of the Polyakov loop as a cusp at 9 = tt/S. Therefore, at 
imaginary chemical potential, the model contains three 
variables; cr, \^\, and the phase (f> as the order parameters 
for the corresponding phase transitions. 
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FIG. 1: The chiral order parameter normalized by as a function of temperature for different values of 6 for the PQM model 
in the mean- field approximation (left panel) and in the FRG approach (right panel). 




FIG. 2: The temperature derivative of the chiral order parameter as a function of temperature for different values of 6 for the 
PQM model in the mean- field approximation (left panel) and in the FRG approach (right panel). 



A. Properties of order parameters in the presence 
of mesonic fluctuations 



By solving the flow equation (|8]) for the quantum po- 
tential we can calculate properties of all three order pa- 
rameters discussed above in the presence of mesonic fluc- 
tuations. 

Figure [T] shows the chiral order parameter normalized 
by = 93 McV as a function of temperature scaled 
by the chiral crossover temperature at = 0, Tpc. To 
identify Tpc the maximum of ~da/dT was used. 

The order parameter is shown in Fig. [1] for various 
values of the imaginary chemical potential. The FRG 
results arc compared to that obtained in the mean-field 
approximation. There is a clear change in the shape of a 
with increasing 6 from to 6 = tt/3. This is the case in 
the mean-field as well as in the FRG approach. However, 
the mesonic fluctuations included in the FRG scheme im- 
ply a strong smearing of the order parameter, making the 
chiral transition smoother. This effect is even more pro- 
nounced in the derivative of the order parameter with 



respect to the temperature. As seen in Fig. [5] the peak 
value of —d{a/ fTr)/dT in the FRG calculations at van- 
ishing chemical potential is almost a half of the mean 
field counterpart. From Fig. [2] one also sees, that in the 
mean-field as well as in the FRG calculations, the peak of 
—d{a/fTr)/dT shifts to higher temperature with increas- 
ing 6. Applying an analytical continuation such a shift 
was used to describe the change of the chiral cross over 
transition temperature with real chemical potential [34 1 . 



In the present calculation, however, one sees two peaks 
structure in —d{a/fTr)/dT. The ones at higher tempera- 
ture are associated with the chiral crossover transition, as 
it is smoothly connected from the unique peak at small 
6. The first peak, which is due to the dcconfinement 
transition, was not previously observed in the mean-field 
calculations in the PNJL model. Its presence in the PQM 
model has to be attributed to the specific implementation 
of the chiral sector. The discontinuity which is specific to 



7r/3, appears atT ~ 1.22 Tpc in the mean field and at 
= 1.39 Tpc in the FRG calculations. This discontinu- 



T 

ity is attributable to the RW endpoint where the phase 
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FIG. 3: The phase of the Polyakov loop (ji as a function of temperature for different values of 9 for the PQM model in the 
mean-field approximation (left panel) and in the FRG approach (right panel). 




FIG. 4: The temperature derivative of modulus of the Polyakov loop as a function of temperature for different values of 6 for 
the PQM model in the mean- field approximation (left panel) and in the FRG approach (right panel). 



of the Polyakov loop shovifs a critical behavior. We note, 
that in the LGT calculations, the RW end point is located 
at T « (1.1 - 1.2)Tpc nil. The obtained model results on 
this temperature strongly depend on the parametrization 
of the Polyakov lo op p otential. For instance, for the loga- 
rithmic potential j32l |. the discontinuity of —d{cr/fTT)/dT 
at = 7r/3 appears at T = l.OSTpc for the FRG approach 
and T = 1.09Tpc in the mean- field approximation. 

Figure 13] shows the as a function of temperature. ^ 
We find in the PQM model with the polynomial potential 
([6]) the second-order transition both in the FRG approach 
and in the mean-field approximation. Thus, the RW end- 
point is a critical point located at 6' = 7r/3 + 27rn/3 and 
at Tptw/Tpc = 1-39 in the FRG approach and at a lower 
temperature T^^/T^c = 1-22 in the mean-field approx- 
imation. Although the order of the transition remains 



^ 9 = -K above Trw is the first-order transition line where two 
minima coexist and Trvv is the bifurcation. For simplicity in the 
figure, we omit the other minimum in 0. 



unchanged by the fiuctuations, in the FRG calculations 
the 4'{T) is a smoother function of temperature. The 
Polyakov loop is treated within the mean-field approx- 
imation also in the FRG approach, thus the observed 
smoothening of (j){T) in Fig. [3] (right) appears owing 
to its coupling to the chiral order parameter which is 
smeared by the quantum fluctuations. 

The effect of the fluctuations on the deconfinement 
transition is more transparent if the derivative of |$| with 
respect to the temperature is considered. The d\^\/dT as 
a function of T is depicted in Fig.|4]in the FRG and mean- 
field calculations. At small 9 this derivative has a unique 
peak, that is attributed to the location of the deconfine- 
ment transition. However, near 6 = 7r/3, where the phase 
(j) also changes rapidly, an additional peak appears. In the 
mean-field calculations one sees two well-separated max- 
ima. At = 7r/3, where the second-order RW endpoint 
appears as the rapid change of seen in Fig.[3J there is a 
discontinuity in d\<^\/dT at T = Trwj both in the mean- 
field and in the FRG calculations. However, the jump 
associated with the RW endpoint is smaller in the FRG 
case. This, as shown in Figs. [T] and [31 is due to mcsonic 
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FIG. 5: Contour plots of the temperature derivative of modulus of the chiral order parameter in the mean-field approximation 
(left panel) and in the FRG approach (right panel). 
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FIG. 6: Contour plots of the temperature derivative of the modulus of the Polyakov loop in the mean-field approximation (left 
panel) and in the FRG approach (right panel). 



fluctuations which smoothcn the critical behavior of the 
chiral order parameter and the phase (p. Consequently, 
as seen in Fig. S] (right), the discontinuity of d|$|/c?T at 
9 ~ tt/S induced by the RW endpoint becomes weaker. 
Nevertheless, in these studies, the deconfinement transi- 
tion, defined by the peak of d\^\/dT, is now separated 
from the RW transition contrary to previous analysis in 
the PNJL model [13, [HI. One should note, however, that 
such separation occurs only in the case of the polynomial 
Polyakov loop potential (|6]), which violates the restric- 
tion on the Polyakov loop target space required by the 
SU{3) Haar measure [51|. If one uses the logarithmic 
potential [32j which exhibits stronger transition, the de- 
confinement critical line connects to the RW endpoint, 
as expected. 



B. Phase diagrams 

The properties of order parameters in the PQM model 
discussed above, can be also represented as the phase di- 



agram in the (T — 9) plane. For illustration, we show in 
Figs.[5]and |6]contour maps of —d{a/fTr)/dT for the chiral 
transition and d\^\/dT for the deconfinement. The effect 
of fluctuations can be clearly seen in the changes of the 
contours. The crossover chiral line, defined by the max- 
imum of —d{a/fTr)/dT, meets the first-order RW transi- 
tion at Tcc/Tpc = 1.47 if the fluctuations are included, 
and at Tce/Tpc = 1.30 in the mcan-fleld approximation. 
However, as discussed above, at imaginary chemical po- 
tential, the transition points for the chiral and decon- 
finement crossover are not determined unambiguously as 
they depend on the particular choice of the Polyakov loop 
potential. In the present calculations, the double peak 
structure can be read off from the contour maps. 

An alternative determination of the chiral critical line 
is based on the location of the minimum in the sigma 
meson mass rricr ■ Figure [7] shows the contour maps for 
TTia, normalized to its vacuum value m'^ . One sees in 
this figure that the double peak structure induced by the 
deconfinement does not show up in the sigma mass. The 
effect of fiuctuations appears as a smearing of the sharp 
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minimum and the stronger curvature. This is evident in 
Fig. [7] when comparing resuhs obtained in the mean-field 
approximation (left) with the FRG approach (right). 

IV. SUMMARY AND CONCLUSIONS 

We have formulated and explored the thermody- 
namics of the Polyakov loop extended quark-meson 
model (PQM) at imaginary chemical potential, includ- 
ing mesonic fluctuations within the functional renormal- 
ization group method (FRG). The flow equations for the 
scale-dependent thermodynamic potential at finite tem- 
perature and at imaginary chemical potential were solved 
in the presence of a background gluonic field. 

We have shown that the non-perturbative fiuctuations 
included in the FRG approach have an important effect 
on the critical properties of the system. Specifically, the 
fluctuations smoothen the chiral phase transition. We 
find that this leads to modification of the curvature of the 
chiral critical line. We also find that the Polyakov loop 
variables (modulus and phase) are also smeared by fluctu- 



ations of the chiral order parameter. We point out that 
the deconflnement crossover transition and the second- 
order Roberge- Weiss endpoint can be separated if the 
polynomial Polyakov loop effective potential is used to 
model the gluon thermodynamics. 

Our results show, that fiuctuations of the meson field 
and non-perturbative effects have to be included when 
studying critical behavior of effective chiral models. Such 
fluctuations can modify substantially the model proper- 
ties and predictions. This is particularly important when 
one attempts to quantify lattice QCD results within ef- 
fective chiral models. 
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